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Abstract 

We have investigated the band- gap pressure coefficients of self- assembled InAs/GaAs quantum 
dots by calculating 17 systems with different quantum dot shape, size, and alloying profile using 
atomistic empirical pseudopotential method within the "strained linear combination of bulk bands" 
approach. Our results confirm the experimentally observed significant reductions of the band gap 
pressure coefficients from the bulk values. We show that the nonlinear pressure coefficients of the 
bulk InAs and GaAs are responsible for these reductions. We also find a rough universal pressure 
coefficient versus band gap relationship which agrees quantitatively with the experimental results. 
We find linear relationships between the percentage of electron wavefunction on the GaAs and the 
quantum dot band gaps and pressure coefficients. These linear relationships can be used to get the 
information of the electron wavefunctions. 
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Self-assembled InAs quantum dots (QDs) grown on lattice-mismatched GaAs(lOO) sub- 
strates have been studied extensively in both experiment and theory in the past 15 years 



due to their potential applications and matured synthesise processes Depending on 

synthesise methods and conditions, the quantum dot can have different size, shape and alloy 
profile. A major task of the research is to study the dependence of the electronic structure 
on the size, shape and alloy profile. The electronic structure includes the electron wave- 
functions and their eigen energies. While there are many experimental ways to probe the 
electron eigen energies and their confinement effects [e.g, photoluminescence (PL) for the 
exciton energyM^he capacitance charging experiment for Coulomb interaction and the single 
particle levels [2]], it is much more difficult to experimentally measure the electronic wave- 
functions. Magnetotunneling spectroscopy Q], low-temperature scanning tunneling spec- 
troscopy 14], and near-field scanning optical microscopes have been used to probe the 
electron wavef unctions, but they are not always successful, and the information about the 
electron wavefunctions remain extremely scarce. Thus any information about the electron 
wavefunctions will be extremely useful. 

One recently popular experimental approach to study the electronic structure of a QD is 
to measure their pressure dependences of the PL energies. While the PL pressure coefficients 
(PC) for both bulk InAs and GaAs are close to 120 meV/GPa, it is found experimentally 
that the PL pressure coefficients for the quantum dots are usually much smaller and they can 
vary significantly from 60 meV/GPa to 100 meV/GPa 0, Q, Is , S 13, Illl depending on the 
samples. While Ma et al. attributed the main reason for the much smaller PC to the built-in 
strain in InAs dots under nonlinear elasticity theory Mintairov et al. emphasized the 
nonuniform In distribution in QD Thus more quantitative analysis and understanding 
are needed here. It is also interesting to find whether the measured PC of a QD can be used 
to infer other properties of the system, e.g., of the electron wavefunctions. 

In this letter, via accurate atomistic calculations for the electron wavefunctions for these 
quantum dots, we show that the nonlinear elasticity and the nonlinear band gap pressure 
dependence are responsible to the reduction of PC. One problem of embedded QD study is 
the lack of reliable experimental information for the QD size and shape. To overcome this, we 
have studied 17 different QD systems covering the possible experimental ranges of QD size 
and shape. What we find, surprisingly, is a universal relationship between the QD exciton 
energy (PL energy) and the pressure coefficients, which can be compared directly with the 
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experimental results. Our calculated PC/exciton energy relationship agrees excellently with 
the experimental measurements. Furthermore, we show that both the QD band gaps and 
their PC correlate linearly with the percentage of the electron wavefunctions on top of the 
GaAs materials. This is independent of the QD size, shape, and alloy profile. As a result, 
these linear relationships and the corresponding PL and PC experiments can be used to get 
the information of the electron wavefunctions. 

We will use the empirical pseudopotential method (EPM) to describe the single 
electron wavefunctions ipi{r) of an InAs quantum dot embedded in a GaAs matrix: 

(-^V^ + ^(r) + VNLmr) = EMr), (1) 

here the total potential V^(r) of the system is a direct sum of the screened atomic empirical 
pseudopotentals z/Q,(r) of the constituent atoms (type a), and Vnl is the nonlocal potential 
describing the spin-orbit interaction. The EPM approach has been used to study InAs/GaAs 
systems extensively, including quantum dots and alloys. Its results agree well with exper- 

various quantum dots in our problem, we need computational 
supercells containing up to one million atoms. The wavefunctions in Eq(l) is expanded by 
planewave basis. In average, each atom will have ~ 50 planewave basis functions. Thus the 
Eq(l) corresponds to a ~ 50 million degree of freedom problem. To solve Eq(l), we have 
used the strained linear combination of bulk band (SLCBB) method In this method, 
the wavefunction ipi{r) is expanded by bulk Bloch states (which is in turn expanded by 
planewaves). Because the bulk Bloch states are good physical basis functions for the quan- 
tum dot states, we can truncate this basis set (down to 10, 000) using physical intuition 
without introducing significant errors. The errors caused by the SLCBB method are around 
10 meV near the band gap compared with the exact solution of Eq(l) ^|. As a result, 
this is a much more accurate method compared to other traditional approaches like the k.p 
method, where a few hundred meV error is possible jl?! . 

To study the pressure effects on the electronic wavefunction, we first need to study the 
. attice relaxation under the pressure. We have used the Keating's valence force field (VFF) 
IS,^^ to described the atomic relaxation. In order to describe accurately the bulk mod- 
ulates and their high order pressure dependence, we have included bond-stretching, bond- 
bending, and bond-angle coupling interactions and high order bond-stretching terms 
Table I lists the VFF bulk modulates and their pressure dependence. They agree well with 
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the experiments. To be able to describe accurately the nonlinear lattice relaxation is im- 
portant because there is a ~ 7.2% lattice mismatch between bulk InAs and GaAs. For a 
quantum dot system, InAs is under compressive stress and GaAs is under tensile stress. 
They will behave differently under additional external pressure because of the nonlinear 
lattice relaxation. 

After the atomic relaxation is described accurately by the VFF model, the pressure depen- 
dence of the bulk band structures for GaAs and InAs is described by the EPM Hamiltonian. 
Here, an explicit local strain dependence of I'aif) is used to describe accurately the defor- 
mation potentials of the band energies Thus the fitting of z/Q,(r) not only provide an 
accurate band structure at zero pressure, it also provides accurate high order pressure de- 
pendence of the band energies. Fig.l shows the calculated band energy pressure dependence 
for bulk InAs and GaAs. The calculated band gap pressure coefficients for InAs and GaAs 
are 117 and 103 meV/GPa respectively, they agree well with the experimental values of 114 
and 106(4) meV/GPa Q. 

We next use the above VFF and EPM Hamiltonians to calculate various embedded quan- 
tum dots under different pressures. A large variety of QD shapes have been reported and 
studied for the InAs/GaAs system by various groups, for example pyramidal quantum dot 
(PQD) with side facets oriented along {101}, {113}, or {105} pi] or truncated pyramidal 



quantum dot (TPQD) [2^|23|. Inside the QD, various In/Ga profiles have been speculated, 
for example an inverted-tri ang le shape In-rich core js^j or a growth direction linearly in- 
creasing In concentration j2J|. To cover the whole spectrum of possible shapes and alloy 
profiles, we have used three sets of QDs: pure pyramidal QDs with {101}, or {113}, or 
{105} facets; pure truncated pyramidal QDs with different height/base ratios; and linearly 
increasing In concentration alloy profile QDs. Besides the shapes and alloy profiles, different 
sizes of the same shape QD are used. In total we have calculated 17 different quantum dots, 
their sizes, shapes and alloy profiles are described in Table II. 

The above described InAs quantum dots are embedded in a pure GaAs matrix. A super- 
cell box is used to contain the quantum dot. A periodic boundary condition is used for the 
supercell box. To remove the possible dot-dot electronic and elastic interactions, sufficient 
GaAs barrier is used. As a result, a supercell can contain upto one million atoms. The 
atomic positions within the supercell are then relaxed by minimizing the strain energy of 
the VFF Hamiltonian. To create a pressure, the overall size of the supercell is changed, and 
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the pressure is calculated from the local GaAs strains away from the quantum dot. After the 
atomic positions are relaxed, the electron and hole eigenstates and eigen energies of Eq(l) 
are solved using the SLCBB method. 

We typically calculate 5 pressure values from to 2 GPa for each quantum dot. Using 
these five points, the band gap of the quantum dot is fitted as Eg{P) = Eg{0) + aiP + a2-P^- 
Then the linear pressure coefficients (PC) of the band gap is read out from ai. In consistent 
with the experiment, we find this PC is in the range of 60-110 meV/GPa, much smaller than 
the bulk InAs and GaAs PC. We then plot all the calculated PC as a function of the QD zero 
pressure exciton energy Eo{0) (which is the band gap minus the electron hole interaction) , the 
result is shown in Fig. 2. Surprisingly, despite all the different shapes and sizes for the 17 QDs 
we studied, we find a rough universal linear relationship between the PC and the exciton 
energy. This provides a convenient way to compare with the experiment, without the need 
to know the QD size and shape which are not available from the experiment. The theory 
and experiment comparison is shown in Fig. 2. The agreement is excellent considering all the 
possible uncertainties involved. We see that, indeed, the QD pressure coefficients are much 
smaller than the bulk values of both InAs and GaAs, and they decrease with the exciton 
energy. 

To understand the variation of the PC, and its dependence on the QD, we can perform a 
simple analysis. We will concentrate on the conduction band minimum (CBM) state since 
most of the band gap pressure coefficient comes from the conduction band [2^. For a simple 
approximation, we can express the energy Ecbm of the CBM eigenstate 'ipcBuij) as a sum 
of an effective mass like potential energy and a kinetic energy Ek-, and the potential energy 
can be approximated by a weighted sum of the local conduction band energy: 



here the E(.{v) is the bulk conduction band energy for the given local strain at r and the 
local constituent material (either GaAs or InAs). Note that, in practice, the space integral of 
Eq(2) is replaced by a sum over the atom WatEc{ai), where the local strain for an atom 
is calculated from the atom's nearest neighbor atomic positions, and Wat denotes the weight 
of |'^CBM(r)P at that atom "at". We have plotted Ecbm as a function of Ylat WatEc{at) in 
Fig. 3 (a). We see that all the calculated QDs fall into a nice curve. The difference between 
this curve and the dashed line (the potential energy line) is the kinetic energy E^. 




(2) 
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Now, we analysis the pressure coefficients of Eqbm using Eq(2). If we ignore the pressure 
dependences of the kinetic energy and the weight function Wat^ we can have an approximated 
relationship: 

E'cBM^ j \i^cBM{r)?K{r)d\ (3) 

here the prime indicate the derivation with pressure. Despite all the approximations, the 
left and right hand side of Eq(3) do form a nice linear relationship, as shown in Fig.3(b). 
The slope of the line in Fig.3(b) is not 1, but 1.25, indicating the right hand side of Eq(3) 
account only for about 80% of the left hand side. This situation can be compared with the 



case of free standing colloidal quantum dots 25[, where the change of PC in a QD can be 



traced back completely from their bulk origin. Our current embedded QD is much more 
complicated due to the internal strain effects between InAs and GaAs, we find such accurate 
analysis is impossible here. 

Despite of not accounting 100% of the left hand side in Eq(3), the physical meaning of 
the right hand side of Eq(3) is clear and useful [especially when it is written as E'^^j^j ^ 
J2at^o-tE'ci^^)]- '^^ quantum dot state is a wavefunction weighted sum of the 

local PC at all the atoms. The E'J^at) depends on the local strain of this atom as illustrated 
in Fig.l. This can be used to understand why the QD PC is in general less than the bulk 
InAs and GaAs results. Because InAs in the QD is always under compressive strain, due to 
the nonlinear PC as shown in Fig.l, the E'^ in the InAs region is significantly smaller than 
its bulk value of 130 meV/GPa. On the other hand, GaAs is under tensile strain, which will 
increase E[. However, because the magnitude of the GaAs strain is in general smaller than 
the InAs strain, and because most wavefunction is localized in the InAs region, the averaged 
PC is then smaller than the bulk InAs and GaAs PCs. Thus we see that the nonlinear bulk 
PC is responsible for the reduction of QD PC compared to bulk values, in consistent with 
the explanation provided by Ma et a/ p. 

Guided by Eqs(2),(3), we now try to find some simple relationships between the experi- 
mentally easily observable quantities (band gap and pressure coefficients) and the wavefunc- 
tion properties. In Eq(3), if we represent E'^{at) by just two values, one for InAs, one for 
GaAs, then E'^j^^j of Eq(3) becomes an linear function of x = J2ateGaAs^a-t/ J2ateaii^at 
(i.e, the percentage of the wavefunction on GaAs). This hypothesis is tested in Fig.4(a), 
where we have plotted the pressure coefficients of the exciton energy (not just the CBM 
energy), so the connection with experiment is more straightforward. We see that E'q and 



X form a very nice straight line. This can be very useful, since a measured E'q value will 
give us the which is a property of the wavefunction that cannot be measured easily by 
other means. The same relationship can be plotted between the exciton energy itself Eq 
and the x, as in Fig. 4(b). They also form a rough linear relationship although with larger 
scatters. The linear relationships in Fig. 4 (a) and (b), in turn, explain why we have a rough 
universal relationship between Eq and Eq in Fig. 2. This is because both Eq and Eq are 
linearly correlated with x. 

In summary, using accurate and reliable empirical pseudopotential methods and the 
SLCBB calculations, we have studied InAs/GaAs quantum dot PL pressure coefficients. 
We investigated 17 different quantum dots covering the ranges of experimental QD size, 
shape and alloy profile. We found a universal PC and exciton energy relationship, which 
agrees excellently with the experimental results. We also find linear relationships between 
the wavefunction percentage on GaAs and the PL pressure coefficient and PL energy. These 
linear relationships can be used to probe the properties of the electron wavef unctions. 
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TABLE I: The VFF bulk modulates and their first and second order pressure coefficients. The 
bulk modulates are in the unit of lOGPa, the dB/dP has a unit 1, and d^B/dP^ is in the unit of 



Property 




GaAs 




InAs 




Fitted 


Expt.'' 


Fitted 




Expt." 


Cii 


12.11 


12.11(4) 


8.328 




8.329 


Cl2 


5.50 


5.48(17) 


4.553 




4.526 


C44 


6.04 


6.04(2) 


3.803 




3.959 


B 


7.70 


7.54 


5.811 




5.794 


dB/dP 


5.01 


4.49 


5.329 




4.787 


d^B/dP^ 


-0.111 




-0.144 







'Reference 121 







TABLE IL 


The 17 calculated quantum 


dots 










Pure InAs pyramidal quant 


um dots 


(PQDs) 








1 


2 


3 


4 


5 


6 


7 




facet 


{101} 


{101} 


{113} 


{113} 


{105} 


{105} 


{105} 




base size(nm) 


6 


11.3 


6 


11.3 


11.3 


15 


20 








Pure InAs truncated 


pyramidal quant 


um dots (TPQDs) 






8 


9 


10 


11 


12 


13 


14 


15 


facet 


{101} 


{101} 


{101} 


{101} 


{101} 


{101} 


{101} 


{101} 


base size(nm) 


6 


6 


6 


11.3 


11.3 


11.3 


11.3 


11.3 


height /base 


2/3 


1/2 


1/4 


2/3 


1/2 


1/4 


1/5 


1/10 



Alloy pyramidal quantum dots 



16 17 
facet {101} {101} 

base size(nm) 11.3 11.3 

alloy profile bottom 40%Ga, tip 0%Ga bottom 50%Ga, tip 0%Ga 
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FIG. 1: The band-edge energies {E(Tqc), and -©(Fst,)) of (a) bulk InAs and (b) GaAs and their 
direct band gap Eg{T8j; — Fgc) under hydrostatic pressure. 
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FIG. 2: The PL pressure coefficient (£"0) versus i?o(0) (PL energy) and comparison with ex- 
periments. The £'0(0) is the zero pressure exciton energy which equals the band gap minus the 
electron hole Coulomb interaction. The experimental results are: Li et al. 11], Ma et al. Ol, Man- 
jon et al. P], and Itskevich et al. 
Williamson et al. 12 1. 
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10 



120 




-4.5-1 , . , . , . 1 -L, . , . , . , r- 

-4.4 -4.3 -4.2 -4.1 80 90 100 110 

s^,W^E^(at) (eV) \^^^'c(^^) (meV/GPa) 



FIG. 3: (a) The Eqbm as a function of X^a^ Wat£'c(ai); (b)The E'^^j^ as a function of 
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FIG. 4: (a)The relationship between E'q{{)) and x (the percentage of the electron state in GaAs); 
(b) The relationship between £'o(0) and x. 
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